The square-lattice spiral magnet Ba2CuGe207 in an in-plane magnetic field. 
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The magnetic structure of Ba2CuGe2 07 is investigated by neutron difFraction in magnetic fields 
applied along several directions in the (a, b) plane of the crystal. In relatively weak fields, H < 0.5 T, 
the propagation vector of the spin-spiral rotates to form a finite angle with the field direction. This 
angle depends on the orientation of H itself. The rotation of the propagation vector is accompanied 
by a re-orientation of the plane of spin rotation in the spiral. The observed behaviour is well described 
by a continuous-limit form of a free energy functional that includes exchange and Dzyaloshinskii- 
Moriya interactions, as well as the Zeeman energy and an empirical anisotropy term. 
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I. INTRODUCTION 

Incommensurate magnetic structures in insulators are 
as a rule caused by a competition between two or more 
magnetic interactions. In the classic description of a spi- 
ral spin structure in Mn02 Yoshimori attributed the spi- 
ral spin arrangement to a competition between two dis- 
tinct antiferromagnetic (AF) Heisenberg exchange inter- 
actions in the crystal.El We shall refer to this kind of spiral 
state as "exchange spiral" . An alternative mechanism is 
realized in a few systems with non-centric crystal struc- 
tures. Isotropic exchange terms in the Hamiltonian com^ 
pete with Dzyaloshinskii-Moriya (D-M) interactionsoEl 
The former are proportional to (Si • S2), i.e., the scalar 
product of interacting spins, and thus favor a coUinear 
state. The D-M term is given by a vector product and 
is usually written as (D • [Si x S2]), where D is the so- 
called Dzyaloshinskii vector that characterizes the ori- 
ented bond between spins Si and S2. D-M interac- 
tions favor a relative angle of 90° between spins. Un- 
like isotropic exchange, D-M coupling is of relativistic 
origin, i.e., is a result of spin-orbit interactions, and ties 
the spin system to the underlying crystal lattice via or- 
bital degrees of freedom. The best known examples of 
real compounds where this second-sgenario ("reLalivis- 
tic spiral") is realized are the MnSicli3 and FeGe.Elcl Re- 
cently we have observed an incommensurate spiral spin 
structure in_B^CuGe207 and explained it within this 
framework.Ell 113 

While the ground states in both models ("exchange 
" and "relativistic") may be similar, under certain con- 
ditions they behave very differently in applied magnetic 
fields. If only isotropic exchange interactions are present, 
the direction of magnetic field is irrelevant to the spin 
structure: the spin space may be freely rotated without 
a change in interaction energy and the spins are not tied 
to any particular direction in the crystal. For "relativis- 
tic" spiral magnets the spin structure in an external field 
H in general may be strongly dependent on the relative 
orientation of H and the Dzyaloshinskii vectors, although 
in several special cases of high-symmetry crystal struc- 
tures this is not so. In FeGe for example, the interac- 
tions of each spin with its neighbors are characterized 
by Dzyaloshinskii vectors pointing in all three equivalent 
orthogonal directions. In this situation applying a very 
small magnetic field always leads to a re-orientation of 
the spin rotation plane normal to the field direction, and, 
consequently, the propagation vector becomes aligned 
along the field. The length of the magnetic propagation 
vector is to a good approximation field-independent. 

In materials of lower crvstaLsymmetry, as fof-|p;am- 
ple in hexagonal CsCuClalla-y and RbMnBrallZllla, or 
tetragonal Ba2CuGe207, the spin rotation plane can not 
be freely re-oriented without affecting the Dzyaloshin- 
skii energy. As a result, the magnetic structure un- 
dergoes drastic changes for certain field geometries, and 
the length of the propagation vector becomes strongly 



field-dependent. An extreme example of this behaviour 
is the commensurate-incommensurate transition that we 
have recently -observed and theoretically analyzed for 
Ba2CuGe207.Eilll3 As a magnetic field is applied along 
the unique c axis of the tetragonal structure, the initially 
uniform sinusoidal spin-spiral is distorted. In a non-zero 
field it may be viewed as a soliton lattice, a regular ar- 
rangement of anti-phase domain wall boundaries separat- 
ing regions of commensurate antiferromagnetic spin-flop 
phase. As the field is increased, the period of the struc- 
ture, given by the distance between solitons, increases 
and diverges at some critical field. He ~ 2.1 T, above 
which a commensurate spin-flop state is observed. 

In this work we address the field dependence of the 
magnetic structure of Ba2CuGe207 in a magnetic field 
H applied in the (a, 6) tetragonal plane. We show that 
the square-lattice spin arrangement in this geometry al- 
lows an almost unhindered re-orientation of the spin 
plane with only a very small change in the period of 
the structure, much like in cubic FeGe. The princi- 
pal difference between FeGe and Ba2CuGe207 is that 
in the former material the propagation vector rotates to- 
wards the direction of applied field. We find that for 
Ba2CuGe207 the propagation vector tends to form a fi- 
nite angle with H, that depends on the orientation of H 
itself. The rotation of the propagation vector with in- 
creasing field is continuous. The observed behaviour is 
found to be in very good agreement with theoretical pre- 
dictions based on a 2-dimensional generalization of the 
simple Ginzburg-Landau energy functional that we have 
previously employed to describe the Dzyaloshinskii tran- 
sition in Ba2CuGe207. 

II. EXPERIMENTAL PROCEDURES 

A single crystal of Ba2CuGe207, the same sample that 
was used in our previous work, was studied in a neutron 
diffraction experiment at the High Flux Isotope Reactor 
at Oak Ridge National Laboratory, on the HB-3 triple- 
axis spectrometer. The sample, roughly 4x4x4 mm'^ 
in volume, was mounted in a small aluminum container 
where it was secured with compressed Al-foil. The 
container was in turn mounted on a precision micro- 
goniometer, that was used to align the (a, b) crystallo- 
graphic plane in the scattering plane of the spectrometer, 
prior to setting up the sample in the cryomagnet. 

The magnetic held was produced by using a split-coil 
horizontal-field superconducting solenoid. The construc- 
tion of the magnet provides two pairs of windows for the 
incident and outgoing neutron beams. The larger pair 
of windows, each window being ±30° wide, is oriented 
perpendicular to the fleld direction. An additional nar- 
rower pair of windows (±15°) is positioned around the 
direction of magnetic fleld. Within the cryomagnet the 
sample could be manually rotated around the vertical 
axis in situ. We estimate the uncertainty in the initial 
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alignment of the sample relative to the field direction to 
be of the order of 1°. After this initial setting the relative 
rotation of the sample in the magnet was done with much 
higher accuracy, « 0.1°, as it could be directly followed 
by monitoring the in-plane Bragg reflections. Measure- 
ments were done for several values of the angle ip be- 
tween the magnetic field and the [110] direction, namely, 
for V =45.0, 50.0, 52.9, 57.5, 60.2, 61.7, 61.9, 71.5 and 
75.2°, the field always being apphed in the (a, 6) plane 
of the crystal. Although the magnet is capable of much 
higher field strengths, stray fields are a serious problem, 
and all the measurements were done in < 2 T. The 
temperature of the sample was controlled by a standard 
He-flow cryostat, allowing us to perform measurements 
in the temperature range 1.8-10 K. As we will discuss in 
detail below, the magnetic structure changes in applied 
fields, but this effect shows a very strong field-hysteresis. 
To avoid this complication all the measurements were 
done using field-cooling. The sample was first warmed 
up to r = 6 K, well aiinve the temperature of magnetic 
ordering T/y = 3.2 K,E] the desired strength and direc- 
tion of the magnetic field were set, and only then was the 
sample brought down to base temperature, T = 1.8 K, 
where most of the measurements were performed. 

At all times the measurements were done in 3-axis 
mode to reduce the signal to background ratio. 60'— 40' — 
40' — 120' collimations were used with pyrolitic graphite 
(PG) (002) reflections for monochromator and analyzer, 
with 14.7 meV neutron energy and a PG filter positioned 
in front of the sample. 

III. EXPERIMENTAL RESULTS 

The crystal and magnetic structures of Ba2CuGe207 
are discussed, in, detail in our previous publications on 
the subject ,c2rO and only the most essential features are 
reviewed here. The magnetcic Cu^^ ions form a square 
lattice in the (a, b) tetragonal plane of the crystal. The 
principal axes of this square lattice, hereafter referred to 
as the X and y axis, are along the [110] and [110] direc- 
tions, respectively. To complete the coordinate system we 
shall choose the z axis to run along [001]. In the ordered 
phase (below T/v = 3.2 K) the spins lie in the (1,-1,0) 
plane and the magnetic propagation vector is (l-f Cj 0) , 
where ( = 0.027. The magnetic structure is a distortion 
of a Neel spin arrangement: a translation along (i, i, 0) 
induces a spin rotation by an angle a — ^ ^ 8.6° (rela- 
tive to an exact antiparallel alignment) in the (1,-1,0) 
plane. Along the [110] direction nearest-neighbor spins 
are perfectly antiparallel. Nearest-neighbor spins from 
adjacent Cu-planes are aligned parallel to each other. 

Obviously, domains with propagation vectors (1 + 
^, ±C, 0) are equivalent and always equally represented in 
a zero-field-cooled sample. This can be seen in Fig. |^(a) 
which shows a contour plot of elastic intensities mea- 
sured at T = 1.8 K on a mesh of points around the 



(1,0,0) position in zero field. We clearly see 4 mag- 
netic reflections symmetrically grouped around the AF 
zone center. The two pairs of peaks, (1 ± C, ±Ci 0) ^^nd 
(1 ± C, =fC, 0), correspond to the two magnetic domains. 
The propagation vector q, measured relative to the AF 
zone center, is strictly along the [110] (or [110]) direc- 
tion. The weak feature around (1, 0, 0) is an artifact that 
is temperature- independent aud was previously identified 
as a multiple scattering peak.li3 Cooling the sample in a 
very weak magnetic field (w 10 mT) always produces a 
single- domain structure. 

When a magnetic field is applied along the [100] direc- 
tion the diffraction pattern begins to change. The prop- 
agation vector q starts to rotate in reciprocal space to- 
wards the direction of applied field [Fig. |l|(b)] , eventually 
becoming perfectly aligned with it [Fig. lJ(c)] . Note that 
all this happens in rather small fields, less than 0.3 T. 
This is to be compared with the field along the c axis. 
He « 2.1 T, that is required to induce the CI transition 
studied previously. Since HcfJ,B ~ |D|, we can conclude 
that for a horizontally applied field the rotation of q oc- 
curs in fields that are an order of magnitude smaller than 
the Dzyaloshinskii energy. 

In the following discussion it is convenient to introduce 
notations for some angles in our experiment (Fig. ||). We 
shall denote by ip the angle between the magnetic field 
H and the [110] direction in the crystal, (p will stand for 
the angle between [110] and q. Using [110] as a reference 
is very convenient since i) = for H = and ii) the 
principal axes of the square-lattice arrangement of the 
magnetic Cu sites in Ba2CuGe207 are along [110] and 
[110], respectively. In these terms, for ip = 7r/4, as H is 
increased, (f> continuously increases from to 7r/4. 

For several directions of the magnetic field we have 
performed careful measurements of (p{H). The results 
are summarized in Fig. ^. Measurements with ip < 61° 
were done on the magnetic satellites around the (1,0,0) 
AF zone-center, and those around (2, 1,0) were used for 
higher values of ip. This was necessary to work around 
the geometrical constraints set by the cryostat windows. 
As mentioned, for ip — n/A the propagation vector starts 
out at = for i7 = and rotates all the way towards 
the field direction {cp = tp = / 4:) , as in the case of FeGe. 
For arbitrary field direction though, while (p levels off 
above H « 0.5 T, it does so at a smaller value, i.e., be- 
fore q reaches the direction of applied field. If we plot 
cp measured at a relatively high field, iJ = 2 T, we find 
that the saturation value for (p is always t:/2 — %p (Fig. 

Thus at high fields the (1,0,0) vector always bisects 
the angle formed by the magnetic field H and the propa- 
gation vector q. Note that in all cases we have ip > tt/A. 
If the magnetic field forms a smaller angle with the [110] 
direction, the magnetic domain with propagation vector 
(1 + OuCfciO), C/iCfc > is destroyed, and the diffraction 
intensity completely shifts over into the OiCfc < domain, 
returning us to the situation when -0 is effectively greater 
than 7r/4. This fact has been verified experimentally. 

An important experimental result is that for all di- 
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rections and of the applied magnetic field the length of 
the magnetic propagation vector |q| is only weakly field- 
dependent. This can be seen in Fig. || that shows the 
trajectory traced by the propagation vector in reciprocal 
space as the magnetic field is increased. For ip = tt/A this 
trajectory is an almost perfect circular arc that starts 
out at (f) — and continues all the way to (f) — 7r/4 
[Fig. 11(a)] . For arbitrary ip we see that \q\ is indeed 
slightly field-dependent and the trajectory has a charac- 
teristic S-shape. This, however, is a weak effect, barely 
detectable with the precision of our experiment, limited 
by the Q-resolution of the spectrometer. 

Having discussed the field-dependence of the propa- 
gation vector, we now turn to that of the spin orienta- 
tion. Typically to solve the magnetic structure one has 
to measure the intensities of several Bragg peak. In the 
present experiment, however, we did not have this luxury. 
The severe geometrical constraints imposed by the nar- 
row cryostat windows eliminate simultaneous access to 
several magnetic peaks. Many of these could be observed 
by rotating the sample within the magnet, however, in 
this case the direction of the magnetic field relative to the 
sample is changed. Only in the case of ip = tt/A, thanks 
to the presence of two orthogonal pairs of windows, could 
we simultaneously observe more than one set of magnetic 
satellites, namely those around (1, 0, 0) and (0, 1, 0). The 
intensities of these reflections are plotted against H in 
Fig. P for -0 = 7r/4. These data are not sufficient to in- 
dependently determine the magnetic structure for each 
value of applied field. We can however assume that for 
all cases, just as for H ^ 0, the spin arrangement is 
a flat magnetic spiral, with the normal to the spin ro- 
tation plane confined in the (a, b) plane of the crystal. 
The assumption is quite reasonable if we consider that in 
our experiments iJ < 2 T and ^mbH < 0.1 meV, so the 
Zeeman energy is at least an order of magnitude smaller 
than the exchange energy (4/5* « 1 meV). A slight con- 
ical distortion of the planar spiral structure is of course 
inevitable in finite magnetic fields, but the tilt towards 
the direction of H will remain very small. If we denote 
by rj the angle between [110] and n, the normal to the 
spin rotation plane {rj ~ tt/2 for H ^ 0), the intensities 
of magnetic satellites around (1, 0, 0) and (0, 1, 0) are de- 
termined by the spin polarization factors in the neutron 
diffraction cross-section: 

^(100) oc 1 + cos^(?7 - 7r/4); 
/(oio) cx 1-f sin2(?7-7r/4). (1) 

From Fig. |^ we see that at high fields /(loo) increases by 
roughly a third of its original value, while /(oio) decreases 
by roughly the same amount. This observation is consis- 
tent with 1] — 7r/4 at high fields, as could be expected: 
a sufficiently strong magnetic field will always align the 
spin rotation plane perpendicular to itself, independent 
of the field direction. For arbitrary tp we therefore expect 
r] = ip in sufficiently high fields. Experimentally, for all 
directions of magnetic field studied, this rule was found 



to be consistent with the observed intensity increase in 
the satellites around (1, 0, 0) or (2, 1, 0) that occurs upon 
increasing the magnetic field from H^QtoH = 2T. 

As will be discussed in detail in the next section, theory 
predicts that for arbitrary values of H one has 

= 7r/2 - 77. (2) 

Since (j){H) is directly measured in our experiments, we 
can plot the expected field dependencies of /(loo) and 
/(oio) using Eqs. (|^) and (||). These are shown in solid 
lines in Fig. ^ and are reasonably consistent with exper- 
imental data. 

We can now summarize all the experimental results 
in a short list of "rules" : i) As the magnetic field is in- 
creased, q rotates in reciprocal space for small H, but 
slows down and eventually stops for H > 0.5 T; ii) |q| is 
only weakly field-dependent; iii) at high fields the [100] 
direction bisects the angle formed by the vectors H and 
q; iv) at high fields n is pointing along the direction of 
applied field; and v) it appears that for arbitrary H the 
[100] direction always bisects the angle formed by the 
vectors q and n. 

IV. THEORY 

We shall now demonstrate that the observed be- 
havior can be well understood using a generalization 
of the approach that we have previously employed to 
quantitatively analyze the Dzyaloshinskii transition in 
Ba2CuGe207. While in our previous studies the direc- 
tion of the propagation vector remained constant, which 
enabled us to use an effectively 1-dimensional model, in 
the present experiment q rotates in the (a, b) plane, and 
we have to extend our expression for the free energy of 
the spin structure to work in 2 dimensions. 

A. The continuous limit 

In general, any almost-antiferromagnetic spin struc- 
ture can be described in terms of the unit vector m(r) 
of local staggered magnetization. This is the continuous- 
limit, when nearest- neighbor spins are almost antiparal- 
lel and, in consequence, the magnetic propagation vector 
is close to that of a Neel structure. For a "relativis- 
tic" spiral, such as the one in Ba2CuGe207, the contin- 
uous approximation is expected to work well when |D|, 
(IbH <^ J , a condition well satisfied in our experiments. 

The simplest form of the free energy that for H = 
gives a planar spin spiral with propagation vector 
along X and the spins rotating in the {x, z) plane, is 
F — Ps/"^ J dxdy{dxm — a/A By x m)^, where ps ~ 
JS^ is the temperature dependent spin stiffness (see 
Refs. [ll| , [T^ ), A = a/\/2 (a being the lattice constant) 
is the nearest neighbor Cu-Cu distance, and a is the 
angle by which spins are rotated at one step of the 
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spiral. For the "relativistic" spiral in Ba2CuGe207, 
a = — arctan |D|/ J « 10°. Of course, on a square lattice 
the spiral can equally well propagate in the direction of 
the symmetrically equivalent y-axis (with spins rotating 
in the {y,z) plane), and a similar term where y is in- 
terchanged with X has to be added to the free energy 
functional. The purely ferromagnetic interactions be- 
tween the Cu-layers in Ba2CuGe207 (Ref. can be ac- 
counted for by including the term f dxdy (Sziii)^, 
where 7 is the spin stiffness anisotropy factor, 7 « 0.03 
for Ba2CuGe207. Finally, the interaction with exter- 
nal magnetic field is given by the Zeeman term Ez = 
-Jdxdy [{x± - ^(H X m)V2 + x\\Hy2] (for de- 
tails see Rcfs. In this last expression x\\ a-nd 
X± are the longitudinal and transverse local staggered 
susceptibilities of the almost AF spin structure, respec- 
tively. For the classical spin model at T = one gets 
Xj_ = ((7/iB)^/(8JA^) and X|| — 0. Combining all of the 
above terms gives us the simplest expression for the free 
energy per Cu-plane, that is consistent with symmetric 
properties of the system and gives the "right answer" for 
H ^0: 

F = J dxdy y ({dxrh - -^e^ x m)^ 
+ {dym - -^Gj, X m)^ + '-fid^m.) 



(X± -X||)(H X rh) 



(3) 



B. Zero magnetic field. 

Let us first consider the system in the absence of a 
magnetic field. We want to find a solution for a general 
direction of the propagation vector in the (x, y) plane. 
For the case when q forms the angle (j) with the x-axis it is 
convenient to change the coordinate system to (x', y', z') 
with x' along the propagation vector, y' perpendicular to 
x' in (x,y) plane, and z' — z. The derivatives and unit 
vectors are changed according to: 

(4) 
(5) 



dx — cos (j) dx' — sin 
dy = cos 4> dy' + sin 
= cos e,xi — sin 

Gyl 



dy, 
dx' 

Gyl 



With these expressions Eq. 
be rewritten as: 



sm (p Gj,/ 
0) after some algebra can 



F 



dxdy 



{{dx'mf ^ {dy,-mf +^{d,'Xhf 



(6) 



2a 

- — (cos 2(f) By, + sin 2(j) ex') • Q^'rh x rh 
2a 

- — (cos 20 ex' + sin 20 Gj,' ) • rh x m - 



+ ^2 (1+^^')) 



(x±-X||)(Hxm)V2-X||i/V2 



For H = 0, since we have selected x' to be the di- 
rection of the propagation vector, dx'Hi ^ 0, while 
9j,'rh = dz'm = 0. It is straightforward to verify that 
the above expression is minimized when i) rh(r) is pe- 
riodic with the period 27r/a and ii) along its vector of 
propagation rh uniformly rotates in a plane that is per- 
pendicular to cos 2(j) ey' + sin 20 g^;' = cos Gy + sin g^,, 
i.e., the "bisection rule" 77 = 7r/2 — formulated in the 
previous section is satisfied. In zero external field all spi- 
ral structures that conform with this bisection rule are 
degenerate, i.e., they have the same free energy. 



C. Non-zero in-plane field 

It is now easy to understand what happens in the case 
H > 0. From all the possible spiral structures, ener- 
getically degenerate at = 0, the system will pick the 
one that takes the most advantage of the Zeeman energy 
— (H X rh)^, namely that which has its spin plane nor- 
mal to the field direction. Independent of the value of 
H (always assuming fisH <ti J), instead of two equiv- 
alent domains seen at iJ = one gets a single domain 
with = tt/2 — -0 and rj = ij). The particular case of 
H||(g^ + ey) is of special interest. In this case the prop- 
agation vector is also directed along (g^, -|- Gj,) and one 
has a "screw-type" spiral with all spins perpendicular to 
the propagation axis. Such a "screw-type" structure is 
realized in MnSi and FeGe for arbitrary direction of q. 

While the above result accounts for the experimentally 
observed behavior for iJ > 0.5 T, we still have to explain 
why the rotations of the spin plane and the propagation 
vector are continuous and in small fields some interme- 
diate structure is realized. Obviously this is due to some 
anisotropy effects that pick the propagation vector along 
[110] (or [110]) and the normal to the spin plane along 
[iTO] (or [110]) for = in the first place. The most 
likely source of anisotropy is spin-orbital interaction. For 
the present discussion however, the actual origin of mag- 
netic anisotropy is of little importance: we can take it 
into account by introducing a phenomenological term 
into Eq. (|^). This is done under the assumption that the 
anisotropy energy Ea, as well as the Zeeman energy Ez^ 
are much smaller than the energy scales of Dzyaloshin- 
skii or exchange interactions. Neither magnetic field, nor 
anisotropy can distort any of the planar spiral structures 
in this limit. Their only effect is to pick the one that gives 
the greatest gain in Zeeman and anisotropy energies. Ea 
can now be written as a function of 77 or 0, which is in 
essence the same thing, since for all the structures we are 
dealing with 77 = 7r/2 — 0_Since at = we know from 
experiment that rj = 7r/2jl3 the anisotropy term must be 
a minimum at this point. It must also comply with the 4- 
fold symmetry of the crystal, i.e., it must be invariant un- 
der ry— >77 + 7r/2. In the most general case it is written as 
Ea = —Ai cos(4?7) — A2 cos(8r7) .... It is also convenient 
to rewrite the Zeeman energy in terms of and t/j: Ez = 



5 



-{X± y XII ) (H X m) 2/2 = ~{X±-X\\)H^ cos^ (77 - ^)/2. 
Minimizing Ez + Ea with respect to rj we obtain: 



^2 _ 8Aisin4?7 + 16A2sin8?7 + 
(X± -X||)sin2(?7- V) 

or, substituting 77 = 7r/2 — 0, 



2 8Aisin40+ 16A2sinJ 



(Xi -X||)sin2((/) + V) 



(7) 



(8) 



This expression gives us the direction of the propagation 
vector and the orientation of the spin rotation plane for 
arbitrary H . 



D. Comparison with experiment 

We can now verify that the analytical results obtained 
in the previous section are consistent with our experi- 
ments on Ba2CuGe207. To begin with, all the approx- 
imations that were made in the above calculations are 
justified. Indeed, the anisotropy energy in Ba2CuGe207 
corresponds to magnetic fields of several tenths of a Tesla, 
fields in which the propagation vector actively rotates. 
The energy of Dzyaloshinskii interactions, as previously 
mentioned, corresponds to fields of the order of 2 T, while 
the exchange energy is roughly 10 times as large. 

Equation (||) can be directly used to fit the experi- 
mentally measured 4>{H). The solid lines in Fig. ^ are 
the result of a global fit of this expression to all our data 
collected at different V'- The only adjustable parameter 
was the anisotropy coefficient Ai = Ai/{xi_ — X||)i while 
we assumed A2 — ~ . . . = Q. The refined value for 
ii is 1.95(0.13) • 10^3 The accuracy of the single- 
parameter fit is quite remarkable. 

It is interesting to note that our theory predicts qual- 
itatively different behaviors for the cases oi ij} = t^I^ 
and '0 > 7r/4, which is particularly easy to understand 
in the case A2 = A^ ~ . . . = Q. For = ''^1^ 
Eq. (|) turns into iJ^ = SAi sin4(?!)/(x± - X||)sin2?7 = 
I6A1 sin 2(/)/(xx — X||)- We see that (j){H) changes con- 
tinuously for ij2 < 16yli/(x± — X||)j has a kink at this 
value and remains constant (equal to 7r/4) above this 
threshold. On the other hand, for tp > tt/A the ap- 
proaches 7r/2 — if) asymptotically in high fields, (/)(i?) is a 
smooth function and there is no threshold field. Precisely 
this kind of behavior is observed in experiment (Fig. ^ . 

The limited data that we have for the orientation of the 
spin plane are also totally consistent with theory. The 
only effect that our theoretical model fails to account for 
is the observed slight variation of |q|. This phenomenon 
may be a result of corrections to the continuous-limit ap- 
proximation that we have ignored. Alternatively, it may 
be a purely quantum effect, related to the small value 
(5 = 1/2) of spins involved and the quasi-2-dimensional 
nature of the system. 



V. CONCLUDING REMARKS 

Again we emphasize the key difference between FeGe 
and MnSi on one side, and Ba2CuGe207 on the other. 
In the former two systems the Dzyaloshinskii vector is 
pointing along the bond between interacting spins, while 
in Ba2CuGe207 D is orthogonal to this bond. As a re- 
sult, in FeGe the normal to the spin plane fi, and the 
propagation vector q are always coUinear. In contrast, in 
Ba2CuGe207 this colinearity is replaced by the "bisec- 
tion rule" : q and n form equal angles with the a axis of 
the crystal, but are not colinear. 

In conclusion, we believe that the present work gives 
a fairly complete picture of the rather exotic static 
magnetic properties of the square-lattice Dzyaloshinskii- 
Moriya antiferromagnet Ba2CuGe207. Much remains to 
be learned in the study of magnetic critical behavior, as 
well as the dynamical properties of the soliton lattice, 
realized in finite fields applied along the c-axis. 
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FIG. 1. Contour plots of elastic neutron scattering inten- 
sity measured in Ba2CuGe207 around the (1,0,0) antiferro- 
magnetic zonc-ccntcr for three values of magnetic field ap- 
plied along the [100] direction. Plot (a) was measured for a 
zero-ficld-cooled sample, and field-cooling was used for plots 
(b) and (c). 

FIG. 2. Experimental geometry: The magnetic field H is 
applied in the (a, b) crystallographic plane at an angle tp to 
the [110] direction. The propagation vector q and the normal 
to the plane of spin rotations form angles (j) and rj with the 
[110] direction, respectively. 

FIG. 3. Ficld-dcpcridcncc of the direction of magnetic 
propagation vector measured in Ba2CuGe2 07 for several di- 
rections of magnetic field applied in the (a, b) crystallographic 
plane. The solid lines represent a theoretical fit to the data, 
described in the text. 

FIG. 4. <j), the direction of magnetic propagation vector, 
measured in Ba2CuGe2 07 as a function of ip, the direction 
of the external magnetic field for H = 2 T. The dashed line 
shows the theoretical result. 

FIG. 5. Measured trajectories traced by the magnetic prop- 
agation vector in Ba2CuGe207 as a magnetic field is applied 
in the (a, b) plane of the crystal. The solid lines are guides 
for the eye. 

FIG. 6. Measured intensities of magnetic satellites around 

the (1,0,0) and (0,1,0) antiferromagnctic zone-centers in 
Ba2CuGe2 07 plotted as a function of magnetic field applied 
along the a-axis. The solid line are the theoretically predicted 
dependences, as described in the text. 
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